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X. . 

We study the spectrum of the Salpeter Hamiltonian H = P^/mF^rY* 2 +V (r) , where V(r) is 
an attractive central potential in three dimensions. If V(r) is a convex transformation of the 
Coulomb potential — 1/r and a concave transformation of the harmonic-oscillator potential 
r 2 , then upper and lower bounds on the discrete eigenvalues of H can be constructed, which 
may all be expressed in the form 



E = min 

r>0 



for suitable values of P here provided. At the critical point r = f the relative growth to the 
Coulomb potential h(r) = —1/r must be bounded by dV/dh < 2/3 /it. 
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I Introduction 

We study the (semirelativistic) "spinless-Salpeter" Hamiltonian 



H = (3^m 2 + p 2 + V{r) , /3 > 



in which V(r) is a central potential in three spatial dimensions. The eigenvalue equation of 
this operator is called the "spinless Salpeter equation." This equation of motion arises as a 
well-defined standard approximation to the Bethe-Salpeter formalism |T| for the description 
of bound states within a (relativistic) quantum field theory and is arrived at by the following 
simplifying steps: 

1. Eliminate all timelike variables by assuming the Bethe-Salpeter kernel that describes 
the interactions between the bound-state constituents to be static, i.e., instantaneous; 
the result of this reduction step is called the "instantaneous Bethe-Salpeter equation" 
or the "Salpeter equation" M. 

2. Neglect the spin of the bound-state constituents, assume the Bethe-Salpeter kernel to 
be of convolution type (as is frequently the case), and consider merely positive-energy 
solutions ip, in order to arrive at the so-called "spinless Salpeter equation" Hip = Eip, 
with a Hamiltonian H of the form (|I]). (For two particles, this form of the Hamiltonian 
H holds only for equal masses m of the bound-state constituents.) 

(For a more detailed account of the reduction of the Bethe-Salpeter equation to the spinless 
Salpeter equation, consult, e.g., the introductory sections of Refs. 0, £§.) This wave equation 
describes the bound states of spin-zero particles (scalar bosons) as well as the spin-averaged 
spectra of the bound states of fermions. 

In this paper we consider potentials which are at the same time convex transformations 
V(r) = g(h(r)) of the Coulomb potential h(r) = — 1/r and concave transformations of the 
harmonic-oscillator potential h(r) = r 2 . The reason for this is that spectral information is 
known for these two "basis" potentials h(r). Thus the class of potentials is those V(r) that 
have a dual representation 



V(r) = gV (-i) = s (2 V) 



in which gQ> is convex > 0) and g^ is concave {g^" < 0). An example of a potential 

in this class is 

V(r) = h C2 In r + c^r + c^r 2 , (2) 

r 

where the coefficients {q} are not negative and are not all zero. Thus tangent lines to the 
transformation function g(h) are of the form ah + b and are either Coulomb potentials lying 
below V, or harmonic-oscillator potentials lying above V. This geometrical idea is the basis 
for our approach to the spectral problem posed by H. We recall the application of this idea 



to the (nonrelativistic) Schrodinger problem in Sec. |HT| . The general envelope formalism for 
the derivation of upper and lower bounds on the eigenvalues of the semirelativistic Salpeter 
Hamiltonian H of Eq. ([!]) is established in Sec. [TV. 



It is fundamental to our method that we first know something about the spectrum of H 
in those cases where V(r) is one of the basis potentials, i.e., the Coulomb and the harmonic 
oscillator. These two spectra are discussed in Sec. [TJ below. In Sec. [V] we look at the example 
of the Coulomb-plus-linear potential. 
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II The Coulomb and harmonic-oscillator potentials 
II. 1 Scaling behaviour 

Since the two basis potentials are both pure powers, it is helpful first to determine what can 
be learnt about the corresponding eigenvalues by the use of standard scaling arguments. By 
employing a wave function 4>(cr) depending on a scale variable c > 0, we find the following 
scaling rule for the eigenvalues corresponding to attractive pure power potentials v sgn(q)r q . 
The Hamiltonian 

H = j3\Jm? + p 2 + v sgn(g)r 9 
has the (energy) eigenvalues E(v,f3,m), where 

The scaling behaviour described by the above formula allows us to consider the one-particle, 
unit-mass special case m = (3 = I initially, that is to say, to work w.l.o.g. with the operator 

H — yl + p 2 + t> sgn(g)r 9 . 



1 1. 2 Coulomb potential 

In the case of the Coulomb potential V(r) = — v/r it is well known || that the Hamiltonian 
H has a Friedrichs extension provided the coupling constant v is not too large. Specifically, 
it is necessary in this case that v is smaller than a critical value v c of the coupling constant: 

2 

v < v c = - . 

7T 

With this restriction, a lower bound to the bottom of the spectrum is provided by Herbst's 
formula 

E >y/l- (at;) 2 , a=|. (3) 

By comparing the spinless Salpeter problem to the corresponding Klein-Gordon equation, 
Martin and Roy || have shown that if the coupling constant is further restricted by v < ~, 
then an improved lower bound is provided by the expression 




■ ■ VI - 4:V 2 1 

£o>\/ — 2 , v<-. (4) 

It turns out that our lower-bound theory has a simpler form when the Coulomb eigenvalue 
formula has the form of Eq. ([3]) rather than that of Eq. @. For this reason, we have derived 
from Eq. (Q) a family of Coulomb bounds which, by rather elementary methods, is found to 
read 



E > y/l - (H 2 , «<^-2<^. (5) 

All these (lower) bounds are slightly weaker than the Martin-Roy bound (fD but above the 
Herbst bound (|3|). We note that these functions of the coupling constant v are all monotone 
and concave. 
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II. 3 Harmonic-oscillator potential 

In the case of the harmonic-oscillator potential, i.e., V(r) = vr 2 , much more is known |7|, ||. 
In momentum-space representation the operator p becomes a c-variable and thus, from the 
spectral point of view, the Hamiltonian H = y/1 + p 2 + vr 2 is equivalent to the Schrodinger 
operator 

H = -vA + Vl +r 2 . (6) 

Since the potential V(r) increases without bound, the spectrum of H is entirely discrete ||. 
We denote its eigenvalues by £ n i(v), where n = 1, 2, 3, . . . "counts" the radial states in each 
angular-momentum subspace labelled by i = 0, 1, 2, ... . Below we shall either approximate 
the eigenvalues £(v) analytically or assume them to be known numerically. The eigenvalues 
£ (v) of such Schrodinger operators are concave functions of the coupling constant v JT0|, ITT]. 



II. 4 The spectral comparison theorem 

For the class of interaction potentials given by (0) with the coefficient of the Coulombic term 
not too large, that is, for all potentials which satisfy the constraint lim r ^ r 2 V'(r) < 2/3/tt, 
the semirelativistic Salpeter Hamiltonian H is bounded below and is essentially self-adjoint 
|J. Consequently, the discrete spectrum of H is characterized variationally || and it follows 
immediately from this that, if we compare two such Hamiltonians H having the potentials 
V^'(r) and V^ 2 '(r), respectively, and we know that V^'(r) < V^ 2 '(r), then we may conclude 
that the corresponding discrete eigenvalues E ni satisfy the inequalities E$ < E^ e . We shall 
refer to this fundamental result as the "spectral comparison theorem." In the more common 
case of nonrelativistic dynamics, i.e., for a (nonrelativistic) kinetic term of the form /3p 2 /2m 
in the Hamiltonian H, a constraint similar to the above would hold for the coefficient of a 
possible additional (attractive) — 1/r 2 term in the potential V(r). 



Ill General envelope theory of Schrodinger operators 



In nonrelativistic envelope theory |TT, [T2|, |T3|], if the potential V is a smooth transformation 
V(r) = g(sgn(q)r q ) of the power-law potential sgn(q)r q , the eigenvalues of H = —A + V(r) 
are approximated by 



E 



nl 



mm 

r>0 



+ V(r) 



(7) 



The numbers P n e(q) can be derived from the eigenvalues of — A+sgn(g)r 9 [^3 ]. If g is convex 
(g" > 0), Eq. (0) yields lower bounds; if g is concave (g" < 0), the results are upper bounds. 

As an immediate application we consider the (nonrelativistic) Schrodinger Hamiltonian 
@ for the Salpeter harmonic-oscillator problem (|J). Here we have H = —vA + V(r), with 
V(r) = (3^/m 2 + r 2 ; hence, the potential is a convex transformation of a linear potential and 
a concave transformation of a harmonic-oscillator potential. We conclude therefore from (|7|): 

1 



mm 

r>0 



Pnl\ 



m 2 + r 1 



< £ne{v) < min 

r>0 



p2 



(8) 



the numbers P n ^(l) are given in Table 1 of Ref. [TJ, and P nl {2) = 2n+i-\. Interestingly the 



upper and lower bounds (§) are equivalent to the corresponding bounds obtained in Ref. 
however, these earlier specific bounds were not derived as part of a general theory. 



4 



IV General envelope theory of Salpeter Hamiltonians 

Let us now turn to our main topic and consider the spinless-Salpeter Hamiltonian of Eq. (|J), 

H = /3ym 2 + p 2 + V(r) , 

and its eigenvalues E. We shall assume that the potential V(r) is a smooth transformation 
V(r) = g(h(r)) of another potential h(r) and that g has definite convexity so that we obtain 
bounds to the energy eigenvalues E. We suppose that the "basis" potential h(r) generates a 
"tangential" Salpeter problem 

H = (3^m 2 + p 2 + vh(r) , 

for which the eigenvalues e(v), or bounds to them, are known. We shall follow here as closely 
as possible the development in Refs. ITT], O, O] for the corresponding Schrodinger problem. 



We express our results in the form of two theorems. 

1. The approximations we shall use from Sec. ||, regarded as functions of the coupling v, 
and also the (unknown) energy functions e(v) of the "tangential" Salpeter problem are 
all concave: e"(v) < 0. The latter result represents the principal claim of Theorem 1. 

2. In Theorem 2 we begin by using an envelope representation for the potential V(r) and 
then demonstrate that all the spectral formulas that follow possess a certain structure. 

Finally, as an application, we specialize to the case of pure power-law "basis" potentials h(r) 
and, more particularly, to the Coulomb potential and the harmonic-oscillator potential for 
which, at this time, we have spectral information (cf. the discussions in Sees. [II.2| and [II. 3|) . 



IV. 1 Convexity of the energy function 

We begin by proving 

Theorem 1 The function e(v) is concave, that is, e"(v) < 0. 

Proof Suppose the exact eigenvalue and (normalized) eigenvector for the problem posed by 
TL = fty/m 2 + p 2 +vh(r) are e(v) and if>(v,r). By differentiating (ip,Hip) with respect to v 
we find e'(v) = (if), hip). If we now apply if>(v, r) as a trial vector to estimate the energy of the 
operator j3^m' 2 + p 2 + uh(r), in which v has been replaced by u, we obtain an upper bound 
to e(u) which may be written in the form e(u) < e(v) + (u — v)e'(v). This inequality tells us 
that the function e(u) lies beneath its tangents; that is to say, e(u) is concave. □ 



IV. 2 The principal envelope formula 

With the help of Theorem 1 we are able to prove 

Theorem 2 (principal envelope formula) Suppose that the operator fi^rri 2 + p 2 +vh(r) 
has the exact lowest eigenvalue e(v), and suppose that the operator (3\frri 2 + p 2 + g(h(r)) 
has the exact lowest eigenvalue E. Then 

E^£ = min[e(f ) - ve'(v) + g(e'(v))\ . (9) 
If g is concave (that is, g" < 0), then E < S; if g is convex (that is, g" > 0), then E > S. 



5 



Proof All the tangential potentials we shall employ have the form V^(r) = a(t)h(r)+b(t), 
where, as in the Schrodinger case, the coefficients a(t) and b(t) are given by 

a(t) = ^M = 9'(Kt)) , b(t) = V{t) - a(t)h(t) = g(h(t)) - g'{h{t))h{t) , 

and r = t is the point of contact of the potential V(r) and its tangent Vw(r). If, for the sake 
of definiteness, we assume that V = g(h) with g concave (i.e., g" < 0), we obtain a family of 
upper bounds given by 

E < e(t) = e{a{t)) + b(t) . 
The best of these is given by optimizing over t: 

E < e{i) = e(a(t)) + b{t) , 
where t, the value of t which optimizes these bounds, is to be determined as the solution of 

e'(g'(h(t))) = h(t) . 

In the spirit of the Legendre transformation ]T5[ we now consider another problem which has 
the same solution; this second problem is the one that provides us with our basic eigenvalue 
formula. We consider 

£ = min[e(f ) — v e'(v ) + g(e'(v))} , 

which is well defined since e(v ) is concave. The solution has the critical point v = g'{e'(v)). 
If we now apply the correspondence h(t) = e'(v), it follows that the critical point v becomes 
v = g'(h(t)) and the tangential-potential coefficients a and b become 

a{t)=g\e\v)) = v , b(t) = g(e'(v)) - ve'(v) , v = v. (10) 

Meanwhile the original critical (energy) value is given by 

£ {t) = e(a(i)) + b(t) = e(v) - ve'(v) + g(e'(v)) , v = v . 

□ 

From the proof of Theorem 2 it follows immediately that, if the exact energy function e(v) 
corresponding to the basis potential h is not available, then, for g(h) concave, concave upper 
approximations e u (v ) > e(v ) or, for g(h) convex, concave lower approximations e\(v ) < e(v ) 
may be used instead of the exact energy function e(v) in the principal envelope formula (^]). 
Then all the lower tangents will lie even lower and all the upper tangents will lie even higher. 
If g is convex, we obtain a lower bound; if g is concave, we obtain an upper bound; because of 
the concavity of e(v ), this extremum is a minimum in both cases. If we wish to use numerical 
solutions to the "basis" problem (generated by h(r)), or if a completely new energy-bound 
expression becomes available, the principal envelope formula (^) is what would be used first. 

Interestingly, in the formula ([|) the tangential-potential apparatus is no longer evident; 
only the correct convexity is required. As in the Schrodinger case Jll| , once we have the basic 
result, the reformulation in terms of "kinetic potentials" is often useful: the kinetic potential 
h(s) corresponding to some basis potential h(r) is given by the Legendre transformation |15| 

h(s) = e'{y) , s = e{v) — ve'{v) . 

Meanwhile the envelope approximation has the kinetic-potential expression V(s) ~ g(h(s)). 

For both the Coulomb lower bounds (|3]) or (|5]) and the harmonic-oscillator upper bounds 
(^) which we have at present, we may express our general results in a special common form 
which will now be derived. 
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IV. 3 The Coulomb lower bound 



We consider first the Coulomb lower bound in which we assume that the potential V(r) is a 
convex transformation V(r) = g(h(r)) of the Coulomb potential h(r) = —1/r. According to 
Sec. 11.2 , in this case all the "lower" e\(v) have been arranged — with the parameters j3 and 
m returned — in the form 



e x [v\ 



(3m 



\ 



From this it follows by elementary algebra that if we define a new optimization variable r by 
e[(v) = h(r) = —1/r, we have 



e\(v) — ve[(v) = (3\ m 2 + 



p2 J 

r> 1 r 

r z a 

Consequently, the lower bound on the energy eigenvalues E of the spinless Salpeter equation 
becomes 



E > min 

r>0 



/ P 2 

p\m 2 + — + V(r) 



v < f3vp 



'IV. 



Here the boundary value vp of the Coulomb coupling v is given, when simply determined by 
the requirement of boundedness from below of the operator ([!]), by the critical coupling v c , 

2 



v P 



71 



and, when arising from the region of validity of our Coulomb-like family of lower bounds (|5]), 
via P = 1/cr, by 

Vp = py/i - p2 < I ( 12 ) 



{P, vp} pairs may be easily generated from the upper bound on the coupling v in Eq. 
The meaning of the Coulomb-coupling constraint is a(t) < fivp, where a is the coefficient in 
the tangential Coulomb potential given by (|10|). 



IV. 4 The harmonic-oscillator upper bounds 

Next, let us turn to the harmonic-oscillator upper bounds. Our main assumption is here that 
V(r) = g(r 2 ), with g" < 0. In this case the only difficulty is that the basis problem h(r) = r 2 
is equivalent to a Schrodinger problem whose solution £ n i{v) is not known exactly. Following 
the discussion after the proof of Theorem 2, let us call the upper bound provided by Eq. (§) 
e n (v) and let us introduce the shorthand notation P n i(2) = 2n+£— 
following parametric equations for e u (t>): 



P. Then we have the 



P 



e n (v) = v— + pVm 2 + r 2 



pr 4 



p2 

<(«) = -7 ■ 



2P 2 Vm 2 + r 2 

By substituting these expressions into the fundamental envelope formula (^) we obtain the 
following upper bound on all the eigenvalues of the spinless-Salpeter problem with potential 
V(r) = g(r 2 ) and g" < 0: 



E ni < min 

r>0 



/ P 2 

pJm 2 + + V{r) 



P = P nl (2) = 2n 



1 
2 



(13) 
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V The Coulomb-plus- linear (or "funnel") potential 

In order to illustrate the above results by a physically motivated example, let us apply these 
considerations to the Coulomb-plus-linear or (in view of its shape) "funnel" potential 

V(r) = -^ + c 2 r , ci > , c 2 > . 

r 

(This potential provides a reasonable overall description of the strong interactions of quarks 
in hadrons. For the phenomenological description of hadrons in terms of both nonrelativistic 
and semirelativistic potential models, see, e.g., Refs. [IB, |T7|.) By choosing as basis potential 
h(r) = — 1/r, we may write V(r) = g(h(r)) with 

g(h) =Cih-j-, 

which is clearly a convex function of h < 0: g" > 0. Thus the convexity condition is satisfied. 
However, we are not free to choose the coupling constants c\ and c 2 as large as we please. It is 
immediately obvious that, for a particular {P, vp} pair, we must in any case have c\ < (3vp. 
For the full problem the coefficient c 2 of the linear term will also be involved. The coupling v 
we are concerned about is given by fllCf). We have 



v = g (e (v)) = — i = = ci + — = ci + c 2 r 




,2 




From this we obtain, for given values of the parameters m and j3 and for a given {P, vp} pair, 
as a sufficient condition for v < (3vp the "Coulomb coupling constant constraint" on the two 
coupling strengths C\ and c 2 in the funnel potential: 

(14) 

In the case {P = l/y/2, vp = 1/2} and @ — m — 1 this condition reduces to Ci + |c 2 < |. For 
Herbst's lower bound (0), i.e., P = Vp = v c = 2/ir, this constraint clearly yields C\ < (3vp. 
There is no escaping this feature of all energy bounds involving the Coulomb potential: the 
constraint derives from the fundamental observation that the Coulomb coupling v must not 
be too large, so that the (relativistic) kinetic energy is able to counterbalance the Coulomb 
potential in order to maintain the Hamiltonian ([!]) with V(r) = —v/r bounded from below. 

For example, if we seek the largest allowed value of the parameter P by solving Eqs. fll2"D 
and ( |T4] ) together, we find that this largest P is given by 



c 2 sin 4 1 2 



cos 2 t(f3 sin t cos t — c{) 



m , P = sint . (15) 



For the Coulomb-plus-linear potential V(r) = —c\jr + c 2 r under consideration, Fig. [1| 
shows the lower and upper bounds on the lowest energy eigenvalue E of the spinless Salpeter 
equation, given by the envelopes of the lower and upper families of tangential energy curves 
(pT|) and fll3f) . In the case of the lower bound (]XT|) , we have used for each m the best possible 



s 



P(m) provided by fll5|). As m — > 0, the "basis" Coulomb problem H = (3\fm 2 + p 2 —v/r has 
energy e(m) — > 0; thus the Coulomb lower bound for a non-Coulomb problem becomes very 
weak for small values of m. Of course, Eq. ( JI3| ) provides us with rigorous upper bounds for 
every energy level. 

In order to get an idea of the location of the exact energy eigenvalues E, Fig. [I] also shows 
the ground-state energy curve E(m) obtained by the Rayleigh-Ritz variational technique || 



with the Laguerre basis states for the trial space defined in Ref. [18|. Strictly speaking, this 



energy curve represents only an upper bound to the precise eigenvalue E. However, from the 
findings of Ref. [18| the deviations of these Laguerre bounds from the exact eigenvalues may 
be estimated, for the superposition of 25 basis functions used here, to be of the order of 1 %. 



VI Summary and conclusion 

In this analysis we have studied the discrete spectrum of semirelativistic "spinless-Salpeter" 
Hamiltonians H, defined in Eq. (|l|) , by an approach which is based principally on convexity. 
We have at our disposal very definite information concerning, on the one hand, the bottom 
of the spectrum of H for the Coulomb potential, h(r) = — 1/r, and, on the other hand, the 
entire spectrum of H for the harmonic-oscillator potential, h(r) = r 2 . The class of potentials 
that are at the same time a convex transformation of —1/r and a concave transformation of 
r 2 includes, for example, arbitrary linear combinations of Coulomb, logarithmic, linear, and 
harmonic-oscillator terms. The envelope technique applied here takes advantage of the fact 
that all "tangent lines" to the interaction potential V(r) = g(h(r)) in H are potentials of the 
form ah(r) + b, and that, by convexity and the comparison theorem recalled in Subsec. |11.4 , 



the energy eigenvalues corresponding to these "tangent" potentials provide rigorous bounds 
to the unknown exact eigenvalues E of H. If e(v) denotes the energy function — or a suitable 
bound to it — corresponding to the problem posed by a "basis" potential vh(r), where v is a 
positive coupling parameter, the envelopes of upper and lower families of energy curves may 
be found with the help of the "principal envelope formula" 

E ~ min[e(f ) — ve'Ou) + g(e'(v))] . 

Here, a sign of approximate equality is used to indicate that, for a definite convexity of g(h), 
the envelope theory yields lower bounds for convex g(h) and upper bounds for concave g(h). 
With the above principal envelope formula at hand, all new spectral pairs {h(r), e(v)} which 
may become available at some future time can immediately be used to enrich our collection 
of energy bounds. If the basis potential h(r) is a pure power, these bounds can be written as 



r>0 



(3\lm 2 + ^ + V{r) 



where the numbers P n e are obtained from the corresponding underlying basis problems. The 
power of this technique is illustrated, in Sec. [V], by our application to the funnel potential, 
V(r) = —ci/r + C2r. For this problem, we have employed both the semirelativistic Coulomb 
and harmonic-oscillator problems to calculate, respectively, lower and upper bounds on the 
energy eigenvalues of the spinless Salpeter equation. 



9 



We expect that such results would provide bounds on the energy eigenvalues for general 
theoretical discussions, or be used as guides for more tightly focussed analytic or numerical 
studies of the spectra of semirelativistic "spinless-Salpeter" Hamiltonians. 
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E(m) 








m 



Figure 1: Lower bounds (L), according to (|TTD, and upper bounds (U), according to (|13|), on 
the energy eigenvalue E of the ground state [(n, £) = (1, 0)] of the spinless Salpeter equation 
with a Coulomb-plus-linear potential V(r) = — ci/r+C2r, for (3 — 1, c\ — 0.1, and C2 = 0.25. 
The lower bound is given by the general result (|Tl~D with the "best" P(m) provided by ([15]). 
In order to satisfy the Coulomb coupling constraint (|15|) , the mass m must fulfil m > a/5/4. 
For comparison, a (very accurate) Rayleigh-Ritz variational upper bound E is depicted too. 



